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ABSTRACT 


This report treats the progress realized during the first year of a 
two-year study aimed at merging spectral and environmental data to 
facilitate the assessment of stress in sunflowers. The study area under 
analysis is in western and northwestern Minnesota along the Red River 
Valley, which is responsible for 90 percent of the state's production of 
sunflowers . 

Weekly ground observations (acquired in I 98 O and 1981) have been 
analyzed in concert with large scale aerial photography and concurrent 
Landsat data. Using multidate supervised and unsupervised classification 
procedures, it has been found that all crops grown in association with 
sunflowers in the study area are spectrally separable from one another. 
Under conditions of extreme drought, severely stressed plants were 
differentiable from those not severely stressed, but between-crop 
separation was not possible. 

The environmental analysis has involved the use of both detailed 
field surveys from the University of Minnesota Crop Pest Management (CPM) 
program and USOA-SRS yield statistics for Minnesota counties. Initial 
regression analyses to estimate sunflower seed yield showed a sensitivity 
to environmental stress during the flowering and seed development stages. 
One of the most important biological factors related to sunflower product! 
In the Red River Valley area was found to be the extent and severity of 
insect infestations. 
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INTRODUCTION 


This report Is submitted in fulfillment of the reporting requirements 
of Contract NAS 9*16427 end covers the research performed thereunder during 
the period from August 24, 1981 to September 30, 1982. This research Is 
being conducted in support of the Early Warning/Crop Condition Assessment 
Project of AgRISTARS and Is aimed at integrating spectral and environ- 
mental data to facilitate esrly determination of stress. The final pro- 
duct of this two-year effort will be a mathematical model based on the 
physiological response of sunflowers to their environment and on the 
manifestations of this response in the appearance of the crop that strives 
to estimate the degree of stress experienced by the crop. 

This report covers only the first year of work to be performed during 
the planned two-year effort. The objectives during the first year were 
to: 1} analyze existing data sets that contain information on the response 

of sunflowers to their environment, and effects of disease, insects and 
yield losses, 2) continue acquisition of field data to incorporate with 
existing data sets, 3) acquire and analyze aerial photographic and Landsat 
data to document spectrally the field conditions observed during the cur- 
rent year, and 4) prepare for the second year of study wherein the afore- 
mentioned data will be synthesized ond field spectral radiometry will be 
initiated as necessary and appropriate. 

Substantive progress has been made in the pursuit of all of the 
aforementioned objectives during the past year. The effort to date has 
been coordinated mutually between the authors of this report with Lillesand 
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Uading tha photographic and Landsat spactral work and Saalay loading 
tha col taction and analysis of tha anvironmantai data sats. Tha structure 
of this raport raflacts this division of effort and is therefore organiaed 
such that wa treat progress to data In tha spectral analysis area in 
Section I, followed by discussion of the environmental data analysis 
activities in Section II. 

I . SPECTRAL analysis ACTIVITIES 
l.i OBJECTIVE/REFERENCE DATA SOURCES 

Tha western and northwestern districts along Minnesota's Rad River 
Valley are responsible for over 90 percent of the state's sunflower pro* 
duct ion. The primary objective of the spectral component of our first 
year of study has been to assess the spectral separability of sunflowers 
from tha other crops grown in this area. In pursuit of this objective we 
have taken an integrated approach of analyzing available Landsat HSS data, 
supporting aerial photography, and ground reference data covering the I98O 
and 1981 growing seasons. Our emphasis has been on data collected during 
1981 (when this study began) in that data from I98O were less available 
and of lower quality, and I98O was a year of extrema drought. Ground 
reference data for the I98I growing season ware provided from several 
sources, yielding locations oT nearly 400 fields (>80 acres) Including; 
sunflowers, potatoes, small grains, sugar beets, pinto beans, and several 
miscellaneous crops of lesser acreage. The University of Minnesota 
Extension Service Crop Pest Management program (CPM) provided detailed 
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information on the condition of epproxiroetely 150 fields of sunflOMers, 
potatoes, and small grains. Centrol , a farm management cooperative, pro- 
vided us with the location of a number of fields of sugarbeets, pinto 
beans and other miscellaneous crops during I 98 I. The 1980 data set 
included some 50 potato and 50 sunflower fields derived from the CPH 
program. 

The CPH program involves field scouts based out of Crookston, 
Minnesota, in the northwest and Morris, Minnesota, in the west. The 
counties covered from each of these district offices are shown in Figure 
I. (Also shown are the limits of coverage for the three Landsat scenes 
needed to covor the CPM districts). The information on crop condition 
monitored by the CPM program is included In a computerized data base 
maintained at the University of Minnesota and includes weekly information 
on the condition of each field, including the presence and degree of 
infestation of various pests, such as insects, weeds, etc. Location, crop 
history, and management practice information for each field is also in- 
cluded in the data base. In short, this data base has proven to be an 
invaluable aid in our analysis of the Landsat data, particularly when 
analyzed in concert with the aerial photography acquired as part of this 
study. This photo coverage has taken the form of large scale (1:10,000) 
70mm color and color Infrared Images of a subset of 60 of the ground- 
monitored fields. The photographs have provided a bridge between the 
Landsat and ground observations and have aided In our substantiation and 
understanding of anomalous conditions appearing In these data sets. The 
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geographic distribution of the photo coverage is shown in Figure 2. 

Figures 3 and 4 depict a representative sample of tho photography. 

1.2 UNOSAT DATA ANALYSIS PROCEDURES 

The Landsat data samples being used in the study are 240x240 pixel 
segments of full scenes which are displayed on an interactive image analysis 
system. This allows the analyst to view the Landsat data in contrast 
enhanced color, electronically magnify the image to observe the full detail 
in the data, and to outline ground areas over which reference data have 
been collected. Forty image segments or "windows" have been analyzed over 
three Landsat scenes in 1980, at three times during the growing season, 
and 46 windows have been analyzed for two dates in 1981. In the remainder 
of this section of this discussion we will highlight the procedures used 
to analyze the Landsat data. In the ext section we will summarize the 
results obtained by analyzing the 1980 and 1981 data sets, respectively. 

1.2.1 Supervised Analysis Procedure 

The 1980 Landsat scenes were analyzed using only a supervised training 
technique, while both supervised and unsupervised training were employed 
to analyze the I 98 I data. The University of Minnesota Image Processing 
Software (UMIPS) system was used to carry out both analyses. In the 
supervised approach, the operator used ground-determined field locations 
and i-inch to the mile county highway maps to carefully outline polygons 
on the interactive display corresponding as closely as possible to the 
corresponding ground derived information. Since the images were geo- 
metrically corrected, a grid overlay showing mile square segments (Public 
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Land Survey sections) was very useful in locating a desired field in rela* 
tion to easily Identifiable natural and cultural features on the base 
map and display. In delineating the polygon to represent the ground 
conditions in a given field the operator used his Judgement as to the 
exact boundaries of the field and avoided edge pixels to the extent 
possible. If anomalous conditions existed within a field as viewed on 
the interactive system, the operator would create a polygon both with and 
without the anomaly. If aerial photographs were available for a field in 
question they were used to help decide the shape and position of the 
polygon to best represent the around observation. The aerial photographs 
were found to be extremely useful in this regard. 

Using the vertices that defined each polygon, the spectral data for 
each field were extracted from the appropriate Landsat tape. To "clean" 
the data base before further analysis, a hi stogram/ range plot was compared 
for all polygons representing the same crop type. Those polygons that 
were noticably different from the rest within their respective groups were 
investigated in an attempt to create representative crop groupings. Samples 
with fewer than 30 pixels were deleted due to their lack of statistical 
reliability. Fields which exhibited high variability, or mui timodal i ty in 
one or several bands were also deleted If reference Information confirmed 
that they were indeed atypical fields. Once the file editing process was 
completed, fields were grouped and named as to their crop type, condition, 
time of year, etc., and a new set of statistics was then created. Two- 
dimensional scatter plots of each Landsat band vs. all others were then 
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prepared and tha pairwlsa transforroad divarganca vatuas anong all ftalds 
wara computad. 

Divarganca is a maasura of tha statisticai saparabiiity of two spactrai 
data sats, and thus, an indicator of tha probability of arror in discrim- 
inating batwaan tham. Divarganca is computad from 

D,j - itr [{:, - - Zj’)] + itr I(zJ' + sj’){V, - Vj)(V, - Vj)'*'] 

whera D|j is tha divarganco batwaan class i and j; tr raprasants tha traca 
of tha indicatad matrix; I is tha covarianca matrix for tha indlcatad class; 
and V is tha maan vector for ths indicatad class. As computad in tha 
UHIPS software, divarganca is scaled to range batwaan 0 and 2000. Ganaraily, 
if tha divergence value between any given pair of spactrai responses is 
in excess of 1500 the responses are considered separable spectrally. Lower 
divergence values indicate overlapping spactrai rafiactancas and consequent 
"confusion" between cover types. 

1.2.2 Unsupervised Analysis Procedures 

In order to more fully explore the spactrai characteristics of tha 
test areas and simulate an operational crop inventory system, an unsupar- 
vised classification was performed on a subset of the Landsat windows for 
1981 . A UMIPS algorithm which is a variation of tha SEARCH program 
developed by tha NASA Earth Resources Laboratory was used for this purpose. 
This algorithm passes a 3x3 pixel moving window over the data set. Any 
window falling below a user-defined maximum variance threshold is analyzed. 
(This operation reduces the occurrence of "edge effect" spectral classes). 
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Thtse s«ts ar« accuntulated and combinad dOMn to a usar-dafinad maximum 
numbar of spactral clasias, with 18 classas baing found to ba an appro- 
priata numbar with which to work. Onca thasa spactral classas wara dafinad 
the various windows wara classified usli s a scaled distance algorithm. 

(A scaled distance measure Is used because covariance data are not generated 
for the spactral classes in the algorithm). Ground rafaranca data wara 
then used to Identify and combine spactral classas into Information 
classes on the Interactive display. With the aid of the UMIPS software, 
colors were assigned to each spectral class on t;ie display using the 
previous supervised training Information as well as the 70ram CIR aerial 
photography, county highway maps, and Information acquired directly from 
farmers in the area. 

1.3 UNOSAT DATA ANALYSIS RESULTS 

1.3.1 Results from Analysis of 1980 Data 

The dates of analysis for 1980 were June 26, July 23. and September 6. 
(The June 26 Image was not available for the northern scene). Again, 
extremely dry conditions prevailed during 1930, causing crops to vary 
widely in their condition. For the mid-summer scene analysis It was 
therefore necessary to stratify each of the sunflower and potato fields 
into two sub-classes based on their infrared reflectance. The sub-groups 
were created using two major criteria. The first was their appearance 
on the interactive display. Highly Infrared reflective fields appeared 
red or pink, while those of low infrared reflectance appeared cyan. Sub- 
sequent analysis of the variance, range of the digital numbers, transformed 
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divergence values, and bispectral plots of each field demonstrated the 
validity of the separate grouping. Due to the extremely dry conditions 
represented in the data set it is entirely possible that the "sunflower" 
or "potato" fields in the low infrared reflectance groups were fields 
under severe stress and/or were simply bare. With this in mind, we 
adopted the same field grouping for the spring and fall scene ane lyses. 
Although we can show separability between types of fields (stressed/bare 
vs. productive) within both potatoes and sunflowers, there was no dif- 
ferentiability between thes<^ classes by crop type, even on a mu i t ! tempore 1 
basis. 


1.3>2 Results from Analysis of 1^81 Data 

The supervised and unsupervised analysis of the 1981 data has been 
completed for only two dates (July 3 and September 19) and for only the 
two most northerly Landsat scenes shown in Figure 1. Due to the extremely 
cloudy conditions realized during the 7981 growing season, these images 
were the only scenes available for analysis. It goes without saying that 
many of the mid-season comparisons we were expecting to make were not 
obtainable due to the poor Landsat coverage— i .e. , one date of mid-sea.w.'' 
imagery over only two-thirds of our study area. 

We have found for the July 9 image tha\ In general, sunflowers are 
separable from potatoes, wheat, barley, alfalfa, and pinto benas. However, 
the separability of sunflowers from sugarbeets on this single date is 
questionable. It should be noted that the number of fields of alfalfa 
and pinto beans was limited in this study due to Inadequate ground information 
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and the Iom percentage of these crops grown In the study area. Therefore, 
our comparisons of sunflowers with these crops are based on a small 
number of samples which were often located in a limited geographic area. 

in addition to showing separability from several other crops, sun- 
flowers also manifested some interesting spectral subclasses. Sunflower 
fields were grouped into geographic regions and growth stages In an attempt 
to account for this spectral separability among subclasses of sunflowers. 
Since each of the windows extracted from the Landsat tapes were approxi- 
mately 8.5 square miles in dimension and were chosen from all major sun- 
flower growing areas contained within the Landsat scenes, the windows 
were used as strata to compare sunflowers among various geographic areas. 
Using these groupings, sunflowers within a specific window were indeed 
separable from sunflowers in many other geographic areas. As expected, a 
north to south gradient was found to exist within a single scene. Possible 
causes for these gradients include a range of soil types, precipitation 
variations, etc. For a iimitec number of fields the CPM program provided 
the growth stages for sunflowers on or within a few days of July 3 When 
fields that had the same growth stage were grouped together it was found 
that several of the different growth stage groups were separable from one 
another. Whereas July 9 is fairly early In the growing season for sun- 
flowers and there was not complete ground cover by the plants at this 
date. It 's suspected that soil reflectance could be influencing our 
spectral analysis. In short, combining the sunflower fields in the above 
two manners did crecte completely homogeneous groupings. In either 
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grouping there wes general, not complete, separability of subgroups. 
Insufficient data at this time was contained within the CPH data base to 
try to check the condition of the crops within the growth stages studied. 

It Is hoped that when this Information Is entered Into the data base that 
this will help explain some of the within-group stage variability. 

To summarize the results of the July 9 analysis, under "normal" mid- 
season growing conditions sunflowers were found to be separable from small 
grains, alfalfa, potatoes, and pinto beans, but were found to be confused 
with sugarbeets. The latter confusion, however, was not a problem In the 
analysis of the September 19 Image data. 

We have found that for the September 19 Image that, in general, sun- 
flowers are separable from sugarbeets, alfalfa, pinto beans, and corn. 
However, the separability of sunflowers from potatoes at this time is 
questionable. As in the summer, the number of fields of alfalfa, pinto 
beans and corn was limited. However, the divergence values of these crops 
vs. sunflowers indicate that these crops are separable, since virtually 
all of the pairwise divergence values are greater than 1500. It must still 
be kept in mind, however, that the fields of alfalfa, pinto beans, and 
corn are few in number and limited In their geographic distribution. 

The separation of potatoes from sunflowers is much less straight- 
forward in the fall. Around half of the divergence values between potatoes 
and sunflowers fall below 1500, indicating that they are not clearly 
separable. It is difficult to say why they are not more clearly separable, 
but perhaps it is due to their senescence. Both sunflowers and potatoes 
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"dry doMn" at this time of the year, thus they may be losing the character- 
istics that earlier in the summer rendered them separable. In addition, 
many farmers use chemicals to artificially dry down both sunflowers and 
potatoes. This may also explain their lack of separability. We do not 
have the information needed to answer this problem fully. It Is an 
interesting problem, however, since the spraying of crops would change 
the spectra] "signature" very rapidly. (Remember that potatoes and sun- 
flowers were separable in the July 9 data set). 

Small grains were harvested by the September 19 image date so these 
fields were probably bare, contained stubble, or replanted to something 
else. In any case, the field locations that were small grains in the 
July 9 image were again separable from sunflowers in the September 19 
Image regardless of what may have been done to those fields. 

Variability still existed within the sunflower data sets themselves 
for the fall imagery. It was noted that relatively higher infrared 
reflectance was received from plants grown In the northern part of the 
state. Although geographic blocking helped to account for some of the 
variability within sunflowers it again did not account for all of It. 

This is of Interest because this "unaccounted for" variability may be due 
to stress. 

The bottom line from the spectral analysis performed to date then is 
that all crops that are commonly grown In association with sunflowers in 
our study area are spectrally separable from one another when at least 
mid-season and fall images from a "normal" year are utilized. Both the 
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supervised end unsupervised anaiyses supported this conclusion. Further 
work is needed to document the spectral manifestation of stress in sun- 
flowers under normal growing conditions. Under drought conditions it was 
shown that severely stressed plants are differentiable from those not 
severely stressed, but between-crop separation is not possible. 
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2.1 UTEMTURE REVIEW 


Sunflower ( Hellenthus ennus L») production hes rapidly increased in 
the U.S. over the past IS years. Prior to 1970. the acreage planted to 
this crop was but a few thousand, hardly significant when compared to 
other agronomic comnod i t i as . Sparked by rapidly developing markets for 
vegetable oil, sunflower acreage expanded to over four million by I98O. 
Host of this increase In acreage has been in the oilseed varieties, which 
are grown as a source of oil and meal. According to Beard (l), in many 
environments, the sunflower yields more oil per unit area of land than 
any other crop. Seed yields generally vary from 1100 kg/ha to 2200 kg/ha, 
but many exceed 3000 kg/ha under irrigation and top management. About 40 
percent of the seed is oil of high quality. 

The world's primary producers of sunflowers are the USSR, USA, 
Argentina, Australia, and Canada, and the agricultural researchers of 
these countries have contributed most significantly to the development 
of improved hybrids and agronomic practices. Since the early I970's, the 
acreage planted to open-pol lenated varieties has been declining, as plant 
breeders in these countries have developed better hybrids. Today, over 80 
percent of the U.S. sunflower acreage is hybrid seed. Hybrid seed has 
lead to the following major Improvements in sunflower production: (a) 
higher seed yield; (b) greater oil content: (c) more uniform attainment 
of seed maturity within individual plants: (d) disease resistance: (e) 
easier combining at harvest: (f) and tolerance to increasing plant 
populations. All of these factors have contributed to the strong positive 
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tr«nd In production statistics ovar the past dacada. Thay IlkaiMisa confound 
any efforts to modal or predict production of sunflower based on historical 
data. For this reason, a general review of the literature on this crop was 
focused on defining the physiological, ecological, and nMnagement charac* 
terlstics that might best direct a screening of Important physical and 
biological variables which could be used as predictors of seed yield. 

This review will describe research findings in the following 4 categories: 
(l) general agronomic and ecological variables Important to sunflower 
production; (2} methods used to score the rate of growth and development 
of this crop; (3) relationships between sunflower phenology and Growing 
Degree Days (GOD); and (4) the crops response to Irrigation and moisture 
deficits . 

Agronomic and Ecological Considerations 

In the mid northern and southern latitudes, sunflowers are grown where 
moisture and/or length of season characteristics may severely limit the 
production of other crops. According to Robinson (22), sunflowers are 
most commonly In rotation with soybeans and small grains. Though not 
considered a drought resistant crop, sunflowers will often produce 
satisfactorily when other crops are seriously damaged by moisture deficits. 
Lindstrom et al (16) report that characteristics accounting for this response 
include: (1) an extensive and heavily branched taproot system exceeding 2 m 
in depth which can extract greater amounts of soil water than some other 
crops; (2) growth and development characteristics which allow tolerance 
of short periods of drought; and (3) a relatively short critical moisture 
period for seed yieid. 
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In the pest yields have been limited by problems with weeds, diseases, 
Insects, and birds. Both Beard (1) and Cobla and Zimmer (5) report significant 
Improvements In weed control and disease resistance. However, because of the 
morphological characteristics of the floral head, the crop remains highly 
susceptible to yield losses from Insects and birds. In the northern latitudes, 
blackbirds and gold finches can consume significant amounts of seed. Univer* 
si ties recommend reducing bird depredation by avoiding bird-prone areas, 
utilizing various scaring devices, or spraying a chemical bird-repellent. 
Cutworms, sunflower midge, sunflower moth, and seed weevils are major insects 
which account for yield losses in the northern latitudes. Table 1 shows a 
typical sunflower insect monitoring schedule for Minnesota. Tillage, weed 
control, and crop rotations are used to keep infestations at low to moderate 
levels. In recent years, the majority of yield losses have been caused by 
sunflower midge and/or seed weevil. According to McBride and Oseto (18), 
damage to individual fields may range from negligible to total loss. 

Chemical control of these and other Insects is difficult. Frequently, they 
feed on areas of the plant that are not readily exposed to chemical sprays. 
Secondly, chemical sprays are used at the risk of disrupting or killing 
needed poilenators such as bees. Bee colonies are frequently established 
by growers to enhance anthesis of their crop. 

The insect problems sometimes confound an interpretation of planting 
date effects. Insects will generally seek out the most susceptible host 
plants, but with year to year variations of insect phenology, it is difficult 
to predict which planting dates will be associated with minimal insect 
Infestations. 

Archaeological evidence supports the belief that common sunflowers 
are native to both Central and North America, where wild species are adapt- 
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td. Th«y may be planted over a wide range of dates and,because of frost 
risk, in general, higher yields are associated with earlier plantings when 
insects or birds do not significantly affect the crop. Both Unger (31) and 
Beard and Geng (2) have researched the response of sunflower hybrids to 
late planting. Cultivars grown in a wide-range of environments respond 
similarly to late planting by accelerating development toward flowering. 

Yield losses can result when sunflowers are sown within 100 days of normal 
frost dates. A commonly used guideline for planting is when the average soil 
temperature first exceeds kS° F in the spring. 

Cheng (A), Cobia and Zimmer (5), and Robinson (22), and others have 
shown that with adequate moisture sunflowers respond to nitrogen fertilizer. 
Phosphorus is frequently applied along with a starter fertilizer. Typical 
application rates of nitrogen, phosphorus and potassium are 50-100 lbs N, 
15-30 lbs P, and 20-k0 lbs K. Sunflowers remove as much nutrients from the 
soil as most other agronomic crops. Placement of fertilizer is usually in 
bands and not directly with the seed, since sunflower germination is quite 
sensitive to fertilizer salts. No micronutrients are recommended and little 
is known relative to their effects on the crop. 

Plant populations under rainfed conditions generally range from 30,000 
to 60,0C0 plants/hectare. With adequate moisture, Robinson et al (25) noted 
a significant increase in yield from 17,000 to 62,000 plants/hectare. However 
tittle effect on yield occurred between 25,000 and 49,000 plants/hectare. 
Beard (1) reports that plant populations for irrigation have exceeded 85,000 
plants/hectare. 

One of the most unique characteristics of sunflowers is the manner in 
which the crop orients its leaves and floral head toward the sun throughout 
the day. This characteristic is termed heliotropism or phototropi sm. Prior 
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CO flowering, the crop creeks the sun throughout the dey, then returns to 
en eastward orientation durinj the night. Following anthesis (opening of 
the flower) the head faces only eastward. Robinson (23) researched this 
phenomena and found that It may be related to row direction of the crop. 
Host sunflower crops are grown with north-south rows rather than east-west. 

Growth and Development Considerations 

SIddlquI et al (28) were the first to develop a widely accepted growth 
and development scale for sunflowers. They devised a descriptive and 
numerical key for delineating the growth and development sequence of the 
crop. Horphological traits are used to Identify both vegetative and 
reproductive stages, from emergence to physiological maturity. The maturity 
stage Is defined as that time when the back of the head turns yellow and 
the phyllaries (small bracts on backside of the head) turn, brown. Sometimes, 
plants are dried with a desiccant, which helps make the seed (or achene) 
moisture content more uniform at harvest. 

The SIddlquI scale Is used to assess yield reductions due to disease 
and insect problems since these pests have differential effects based on 
the crop stage of growth. A modified version of this scale is used by 
hail insurance adjustors to settle claims. Schneiter and Jones (2f) 
reported that defoliation by hail at the bud or early flower stages 
(3.1 to 4.1 on the SIddlquI scale) can lead to a 20 to 30 percent reduction 
in seed yield if other environmental factors are not limiting. They found 
that hail damage In the vegetative stages did not significantly reduce 
yields. 

Marc and Palmer (17) developed a similar scoring index for sunflower 
development to use in quantitative analysis of the relationship between 
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•nvIroniMntal var{«b1«s and sunflowar ontogany. Soma cuittvars hava thoMn 
a photoparlod rasponsa In thair davalopmant, whila many othars appaar to 
ba day-nautral. 

Coultas (6) raportad on a modi f lad davalopmant Indax for sunf lowers , 
with numerical, verbal , and sketched descriptions of each stage. These 
appaar as Table 2 and Figure 5» respectively. This indax Is used frequently 
by field scouts In Integrated pest management programs. Note that this Index 
describes the ontogeny of a 115*day hybrid and Infers a certain amount of 
stability by the relationship to growing degree units. This will be discussed 
in the next section. 

Relationship to Growing Degree Days (6D0) : 

Robinson et al (20) were some of the first researchers to study the 
relationship between sunflower development and daily temperature accumulations. 
They found a strong relationship between the rate of sunflower development 
and the accumulation of temperature above a base of 7*2^C (k5*^F) . The 
computation of daily Growing Degree Days using this base is shown as: 

GDD • daily mean temperature-base 

where the base equals When this yields a negative quantity, then 

GDD are assumed to be zero for that day. They found that cumulative GDD 
related closely to the rate of development for similar sunflower cultivars. 
Furthermore, they tested this system for latitudes ranging from 31 to 49 
degrees. Each degree of latitude increased the length of the planting to 
ray flowers period by 2 days. However, the GDD for this period varied in 
a narrow range from 829 to 948 from Lincoln, Nebraska to southern Manitoba. 

This indicates some degree of ontogenetic stability in individual cultivars 
when they are grown in wide ranging latitudes. The data suggest that the 
GDD relationship to development rate may be more stable for day-neutral 
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Pitid Ot<wlop « Wi > t of ttw Sun fl onw r Tabl* II. 


•taga 


OggaeipBoii 

1. Emarganca 


From saading until opposlM Navas ara fully axpandad. Each laaf staga is 

SN‘ 

1.1 -1.4 

distinguighad from anotfwr whan tha laaf padola it vigibla ttirougfi tha crown. 

OFF* 

0-14 


GOU' 

0-1B0 


2. Vagatativa 


From tha formation of tha first altamata laaf until laaf formation oaasM. 

SN 

2.1 -2.4 

Each suooatdva altamata laaf aftar tha formation of tha 4th altamaw laaf can 

OFP 

14-60 

ba a substaga of this stage. Total number of laavas can ranga from 26 to 38. 

GDU 

1BM239 


3. Bud 


Terminal bud forms a haad rathar than a duttar of Navas. Plant has raaehad 

$N 

3.1 -3.4 

maximum height and a flower bsxf it distinctly visibN (% to 1 inch buttonl. 

OFP 

60-70 

Inflorascanoa taparatas from Navas and bagint to open. 

GDU 

1240-1510 


4. Anthasis 



SN 

4. 1-4.5 


Early Bloom 

PataN of ray flowers become vNibN as inflorascanoa opens. First anthers 

SN 

4.2 

appear at tha outer edga of tha infloratcanca as anthasis begins in outar 

OFP 

70-60 

circumNranoa of infloratcanca and prograttas to ona4taif Woom. 

GOU 

1511-17B0 


50% Bloom 

Ona-hglf of tha disk flowers are in anthasis. Saad fillirtg begins in 

SN 

4.3 

ouNr florets artd progratsn inward. 

OFP 

80-00 


GOU • 

1781-2040 


Full Bloom 

Anthasis it compNta. Saad filling continues. Seeds at outar edge of 

SN 

4.5 

haad become dark in color. 

OFP • 

90-100 


GOU • 

2041-2280 


5. Saad Oaval 

Bpmant 

Patels of ray flowars wilt and drop from haad. Saad filling continues. 

SN 

5.1 -5.2 

Lower Navas begin to sottesca. Haad becomes irtvertad and is green in color. 

OFP ■ 

100-114 


GOU • 

2281-2545 


Physiological Maturity 

Back of haad it yallow and bracts ara brown. Stems arKl saads ara hard 

SN 

5.3 

and mature. 

OFP • 

115 


GOU • 

2545-2561 



'Growth itfltit Mcordint to tiddtqui. Srown, Mwt Alton r*nont OiioMi Itaportor," t97S. Voi. S9. pp. 7-11). 

*Ooyt from planting lor a fioM to roach eomplation of porticwiar growth tiaga. 

*Appia>ifnata g*owir«g Pagraa unit accumulation for a fioip of lunflowar to roach a gwan ataga of growth. QOU calculatad from a bam 
tamporatuio of 4S* P for a Ilg^V hybrid. 


PIgttfg 

Stagg l.t 
Sttgg 1.2 
Stggg 2-2.4 
(vogatativo) 
Stagg 3.1 
Stagg 3.3 
Stagg 3.3 
Stagg 4.0 
Stagg 4.1 
Stagg 4.1 
Stagg 4.3 
Stagg 4.4 
Stggg 4.5 
Stagg 5.1 
Stagg 5.1 
Stagg 5. 15 
Stagg 5.3 


noaarintktn 

CotyMom gmorgad, patiolaf of fint oppoutg Iggvai not visiWa 
First pair of oppogita laavgi Pgvaiopa J . patielai of laeond pair not visibia 
First aittmata laaf to last aitamatg loaf (dratwing of all laaf stagMl 

Hggd visiblg, % incti “bonon" 

Bud glavatod abova crown taavas 
Croit-saetion 

First antttgrs visiblg on outar odga of infloroscanca 

Early bloom 

Crosa-saction 

50% bloom saad filling in outar florats 
75% bloom 

Full Woom-anthasis eemplatg, saads in outer florats dark in color 
Haad inverts, petals of ray flowars dros, from hggd 
Crou-saction 
S0%mgturg 

Physiological maturity; back of head Turin yellow, brgctt brown 
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cuittvars than thosa showing g. photoparlod rasponsa. Indaad, Goyna at al 
(12) showad that a GDD systatn workad wall in pradicting sunflowar davalop 
ratas In Australia, with tha axcaption of tha flowaring pariod of eultivars 
with strong photoperiod response. 

Robinson (21) in a separata study showad that tha base 7.2®C GOD works 
reasonably wall for different hybrids grown In tha northern latitudes. For 
six common varieties, he found that tha range in maturity as defined by 
total GOO from planting to complete sanascanca was more narrow than that 
reported for corn hybrids, ranging from about 1250 to 1350 GDD. Vary late 
plentings had somewhat lower GDD requirements to reach maturity. GDD were 
shown to be better predictors of development than days from planting for 
al 1 6 eultivars. 

Coultas (6) used Robinson's work as a basis for developing a Fahrenheit 
GDD system (base i»5®F) for typical Minnesota sunflower hybrids. This system 
is used to estimate the develop of sunflowers such that schedules for insect 
scouting programs and irrigations can be better planned. Note in Table 2 
that Minnesota sunflower hybrids typically need approximately 2550®F GDD 
to reach physiological maturity. 

Hammer et al (1^) used the GDD-phenology relationship to develop predictor 
variables for forecasting sunflower yields. They used a moisture stress 
index (based on evapo transpiration) during the GDD predicted flowering 
stage to forecast yields. The model validated very well for one cultivar 
grown in experimental plots. Moisture stress and yield modeling are discussed 
in the next section. 

Response to Irrigation and Moisture Stress 

Levitt (15) defines sunflower as a mesophytic plant with some drought 
tolerance and a water requirement (water used In transpi rat lon/dry matter 
produced) ranging from 250 to 350. This water requirement is somewhat less 
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than that of other mesophytas such as the small grains, but more than that 
of more xerophytic species such as sorghum. Ennen et a I (9) studied the 
comparative >uter use of 8 common oilseed cultivars and other major crops 
grown In the Red River Valley along the Minnesota-North Dakota br.*der. 

Total seasonal water use among cultivars and across locations only ranged 
from 30 to 34 cm (12 to 13.5 inches). Figure 6 shows a comparison of sun- 
flower water use and that of other major crops In the Red River Valley 
based on data from 1973 to 19'7. Note that sunflowers averaged 13 to 14 
inches while alfalfa and corn required an average of 24 to 21 inches 
respectively. 

Lindstrom et al (16) showed that under limited rainfall conditions, 
seasonal water use of sunflowers can be nearly comparable to that of corn, 
however sunflower production has been successful in areas where limited 
moisture has excluded corn production. They explain that the ability of 
the sunflower to produce a crop with limited moisture is not due to a lower 
water use, but to the crop's ability to escape or tolerate drought situations. 
Sunflower appears to be drought tolerant even for the three week interval 
from heading to flower completion, when it exhibits an ability to halt flower 
development during extreme moisture stress and resume when adequate moisture 
is restored. Other crops, such as corn, are susceptible to drought for 
longer periods and to a greater degree during their reproductive stages. 
Indeed, Seeley and Spoden (27) showed that peak daily evapotranspiration 
rates of sunflower, though similar to corn, are not maintained for nearly 
as many days during the reproductive and filling periods. 

The sensitivity to moisture during the flowering phase, though not as 
responsive as corn, is nevertheless a factor of importance in irrigation. 
Sior.it (29) reported that prolonged moisture stress (two weeks or nore) 
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during flowering can reduce seed yields by 35 to 40 percent. For other growth 
stages, yield reductions ranged from negligible to 25 percent, with the 
least effects associated with the vegetative and senescence phases. Using 
this sensitivity in another way, Goyne et ai (13) and Unger (30) reported 
that adequate moisture (via irrigation or rainfall) for maintaining near 
potential transpiration during tne bud and flowering phases boosted yields 
significantly. Cobia and Zinmer (5) found that sunflowers respond as well 
to irrigation as alfalfa and soybeans, but not as well as corn. The yield 
response to irrigation varied little over plant populations ranging froir 
12,000 to 25,000 plants/acre. 

Maximum leaf photosynthetic rates under irrigation at the flowering 
stage range from 38-40 mg C 02 dm hr for sunflower, quite similar to other 
species, but considerably less than Ct^ crops such as corn, which has 
shown maximum rates of 55-60 mg C02dm hr . Boyer (3) suggests that 
recovery of photosynthesis following moisture stress may be more rapid In 
sunflower than many other agronomic crops. 

Water use efficiency (WUE) of sunflower has been researched to a 
limited extent. Unger (30) reports WUE as low as 25 kg ha”^cm”^ under 
moisture limiting conditions. However, Cox and Joliff (7) reported WUE as 
high as 105 kg ha ^cm ' for cultivar 894 when grown in a favorable modified 
marine climate with very low vapor pressure deficits (1.5 KPa) . 

Ennen (10) reported WUE ranging from 37 to 53 kg ha ^cm ^ for sunflower 
yields of 1200 to 1500 kg/ha obtained under rainfed conditions in North 
Dakota. Irrigated yields of 3000 to 3700 kg/ha had WUE up to 76 kg ha”'cm”‘. 
Cheng and Zubriski (4) found that irrigation and increasing nitrogen 
fertilizer applications up to 56 kg/ha raised yields from 2700 kg/ha to 
4200 kg/ha and increased WUE by 25 percent. 
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No significant resaarch or documantatlon was found on the affact of high 
tafflperatura strass on sunfiowar saad ylalds. Hammar at al (14) and Fararas 
at a) (It) reported on methods of relating evapotranspiration (ET) to yields 
and most ET computations consider temperature as a driving variable. 

In summary, the literature showed that the following 7 major characteristics 
are important considerations for sunflower research and modeling efforts. 

1) Sunflower is a native specie of Central and North America. 

2) Sunflower is a C 3 plant and behaves similarly in many ways to other C 3 
species. 

3 ) Hel iotropism, the characteristic of facing the sun until anthesis occurs 
is a unique trait with respect to most other agronomic crops. 

k) Unlike field corn, sunflower hybrids exhibit a more narrow range of 
maturity types. 

5) Sunflowers still benefit from an abundance of pollenator activity such 
as bees. 

6 ) Because they are grown in latitudes with marginal growing seasons, 
sunflowers are sometimes sprayed with desiccants to enhance the 
uniformity in drying of the seed head. 

7) Though progress has been made in both disease and weed control, 
sunflowers are still very susceptible to insect and bird damage. 
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2.2 DATA 


The University of Minnesota Crop Pest Management (CPM) program conducted 
a detailed pest survey over thousands of acres of Minnesota cropland during 
the past 3 years. Sunflowers, being one of the major commodities produced in 
the state, were heavily surveyed in the west central and northwest districts 
of the state, shown in Figure 1. Weekiy survey reports were filed by field 
scouts with the University of Minnesota, where they are stored as part of a 
pest management data base. Sample A shows an example of the weekly field 
survey form used for the 1980 and 1981 crop seasons. 

A subset from the 1980 and 1981 sunflower data base contained a list of 
fieids where end-of season yield samples were taken. Combining the data from 
1980 and 1981 produced a total of I30 yield reports for rields ranging in 
size from 20 to 80 acres. These fields also had continuous weekly survey 
reports from CPM field scouts. This data set was the basis for developing 
environmental variables which could be used in regression analysis with yields. 

Environmental variables were derived from daily weather observations by 
National Weather Service cooperative observers and by researchers at the 
University of Minnesota Agricultural Experiment Stations at Crookston and 
Morris. These observations contained records of daily maximum temperature, 
minimum temperature, and precipitation. Environmental data from sixteen 
different observers were used to describe the growing conditions for sunflowers 
In the 130 CPM fields. This represented a 13 county area located along the 
Minnesota border with North and South Dakota. 

For the original 130 fields, a screening procedure eliminated 29 of 
them. These fields were deleted from the data set for various reasons 
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SAMPLE A 


Ac^icuttursl CKMntion Sarvict 
UniMrtitv of Mlnnatota - U.S. Oapartmtflt Agrieuliur* 
tnitituM of Atricultwra. PorMtry and Homa Eeenomict 
St. Paul, MiiwaaotaSBIOS 


Field Survey Form - SUNFLOWERS 


I. FIELD INFORMATION 

Coopvrator't ntmt 

Cooperator's field I.D. . 

Scout name & 1.0 

Date 

II. CROP DESCRIPTION 

NW SW C SE NE 


e e e • 

e e e e 
a a a e 
a a a a 


1 . Plant height (average in inches) 

2. Plant population (1/100 A) 

3. Planting titwa 

4. Growth stage 


III. INSECT OBSERVATIONS (Total insectt/20 plants at each location or % defoliation) 

NW SW C SE NE Avg. Total 

1. Stem weevil 

2. Sunflower beetle (%defol.) 

3. Sunflower moth - 

4. Painted lady (% defol.) 

5. Sunflower midge . — 

6. Seed weevil 

7. Banded sunflower moth . ■ 

8. Grasshopper - - 

9. Other - — 

Comments — 


IV. SUNFLOWER DISEASES - Severity (0-9 scale) & % Incidence 

NW SW C SE NE 


Distance 

Walked 


No. of Diseased % 

Plants Inadanca 


1 . Downy mildew 

2. Sdarotinia 

3. Phoma 

4. Fusarium 

5. Rhizopus 

6. Verticillium 

7. Other 

Comments 


V. WEEDS 

Common Name 


Code 


Stage 


Severity 


Comments 

VI. CROP INJURY 

□ Yes GNo Location in Field: GnW DswOcGseGnE 

VII. COMMENTS 
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Including Irrigation, quastlonabla ylald astimatas, Incomplata Moakly survay 
reports, disaasa or insact Infastatlons not accountad for In tha survays, or 
lack of raprasantative anvironmantal data from a location In closa proximity 
to tha flald sita. 

Tha data from the remaining 101 fields, along with computed anvironmantal 
variables are shown in Appendix I. Tha data are separated by blank colixnns. 
Thera are a total of 16 records In each row of data. Table Ml shows tha list 
of variables used in tha original regression analysis on tha CPN sunflower 
data base. All data were derived from the field survey reports and local 
weather observations. Growth stage data ware collected using a modified 
Siddiqui type scoring system, where bud, anthesis (bloom), and seed devel- 
opment phases were denoted as stages 3.1, 4>1» ond 5.1, respectively. Potential 
evapotranspi ration (PET) estimates were calculated by the USDA-SCS Blaney- 
Griddle Method, which uses local daily temperature and day length data along 
with empirically derived crop coefficients for sunflower growth stages. Field 
capacity coefficients (FCC) are based on soil texture data and represent 
relative water holding capacities of the soils. The insect codes pertain only 
to infestations of sunflower seed weevil and midge, which were found to be 
the primary causes of yield reductions in the I98O and I98I crop seasons. 

Insect counts and thresholds are based on number of larvae found on 20 plants 
In 5 different samples from each field (See Sample A). 

The CPM data base provided the opportunity to study the variations in 
yields from individual fields which represented a range of soils and climate, 
but only spanned 2 crop seasons. In order to look at yield response from a 
different perspective, USDA-SRS county yield data from 4 Minnesota counties 
for the period 1971 to 1980 were used in a separate regression analysis. 

1974 was excluded due to severe drought and crop failure. These data are from 
Marshal), Polk, Norman, and Traverse counties where the acreage planted to 
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Table Ml 

Variables used in regression analysis on the CPH sunflower data base for 
1980 and 1981. 


YD ■ Yield in Ibs/acre. 

YR ■ Crop season. 

FCC ■ Field capacity coefficient or 

relative water holding capacity 
of the soil. 

PD * Planting date (Julian day). 

ft 

PP * Preseason precipitation (inches) 
from April 1 to May 31 • 

AD > Anthesis date (Julian day). 

FP « Precipitation (inches) from 10 days 

before to 10 days after flowering. 

PET3 “ Potential evapotranspi ration (inches) 
from 5 days before to 5 days after the 
bud stage. 

P3 • Precipitation (inches) from 5 days before 
to 5 days after the bud stage. 

PETA ■ Potential evapotranspi ration (inches) from 
5 days before to 5 days after anthesis. 

PA ■ Precipitation (inches) from 5 days before 
to 5 days after anthesis. 

PET5 • Potential evapotranspi ration (inches) from 

S days before to S days after seed development. 

P5 • Precipitation (inches) from 5 days before 
to 5 days after seed deveiopment. 

FPP ■ Final plant population (thousands/acre). 

IC " Insect code for seed weevil and midge, 

indicates relative degree of infestation. 
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sunflowers is the highest In the stete. Environmental date were collected 
from local National Weather Service cooperative observers and averaged for 
each county. These data appear In Appendix II, where they are defined. Since 
growth stage data were not available, with the exception of USDA-Crop Reporting 
Service estimates of planting dates. It was not possible to compute environ^ 
mental variables similar to those used In the CPH data base. Several other 
detailed field data were not available, and therefore the predictor variable 
set was quite small. However, because of the longer time series of yields, 
both time trend and planting date effects could be studied. Table IV lists 
the variable names used in regression analysis of this data set. Note that 
the precipitation and temperature variables are oriented to an estimated 
date of anthesis, which was computed based on GOD of base (after Coultas). 
The temperature stress variable is used- In this analysis as a substitute for 
the PET variable computed for the CPH data base. For the county data sets, 
it was not possible to accurately use the USDA Blaney-Criddle Method to 
estimate PET. 

A discussion of the regression analyses on these two data sets follows. 
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Table IV 


Variables used In regression analysis on the county sunflower yield data 
from 1971 to 1980, excluding 1974 . 


CAYD ■ County average yield In Ibs/acre. 

TTV ■ Time-trend variable, where. 

1 - 1971 and 9 - 1980. 

CAPO > County average planting date, when 

50 percent of the acreage was planted 
(noted as Julian day). 

CAPF ■ Average county precipitation from 10 

days before estimated date of flowering 
to 10 days after (Inches). 

CATSF ■ County average of number of days when 
the daily maximum temperature equaled 
or exceeded 90°F during the flowering 
period. 
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2.3 MODELS 


Multiple regression analyses were dona on both the CPM and county data 
sets. All variables and variable combinations were screened based on simple 
correlation coefficient, contribution to R , t*test, and F-test for each 
model. Only the best regression models are discussed In this section. Further 
analyses could be done by the reader using the data presented in Appendices 
I and II along with any of a number of available statistical analysis software 
packages. 

Regression Analysis on the CPM Date Set 

Tables V and VI summarize the statistics from the best fit model for 
yield using the CPM data base. The 1C variable was by far the most significant 
predictor, accounting for 62 percent of the yield variation. Of the lOt cases, 

36 showed some yield damage from seed weevil, midge or both. This represents 
36 percent of the sunflower fields scouted and may or may not be indicative of 
the important role insects play in northern latitude sunflower production. 

Both 1980 and 198I are described by entomologists as relatively high activity 
seasons for sunflower insects, and may not be most representative of average 
conditions based on more years of data. The PETS variable was the second 
best predictor, indicating that the transpiration demand on the crop during 
the seed development stage, at least, partially determines yield. This would 
agree somewhat with findings reported in the literature. The only other variable 
significant at P ■ .01 level was PP, indicating the importance of early season 
stored soil moisture. 

2 

The other two variaMes, of lesser importance, were FP and FP (FP squared), 
both significant at P ■ .05 level. This relationship agrees very well with 
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TabU V 


A sumniry of simpto statistics from tha bast ragrasslon modal for sunflOMar 
ylald basad on tha CPH data basa. 


Varlabla 

(units) 

Mean 

Value 

Standard 

Deviation 

Coefficient 
of Variation 

Maximum 

Value 

Hinimum 

Value 

YD Obs/acra) 

1ii3A 

311 

.217 

1800 

600 

PP (Inchas) 

5. 84 

1.53 

.262 

8.98 

3.89 

FP (Inchas) 

3.27 

1.53 

.468 

7.19 

1.1*7 

PETS (Inches) 

1.94 

.20 

.104 

2.39 

1.44 

1C (coded) 

0.42 

-60 

1.453 

2.0 

0 

FP^ (Inches) 

13.01 

13.02 

1.000 

51.7 

2.16 


Tabla VI 


A summary of 
the CPH data 

regression statistics 
base. 

for the 

best sunflower yield 

model from 

Variable 

Regression t 

Coefficient 

-value 

Contribution 
to R^ 

Stepwise 

r 2 

Intercept 

1822.37 

11.58 

- 

- 

PP 

20.07 

2 . 78 ** 

.07 

.07 

FP 

80.39 

2.00** 

.07 

.14 

PETS 

-232.71 

- 4 . 23 ** 

.13 

.27 

1C 

-440.62 

23.14** 

.62 

.89 

FP^ 

- 10.37 

- 2 . 19 * 

.01 

.90 


*SlgnlfIcant at P ■ .05 Of ■ 95 

**Slgnlf leant at P ■ .01 
RMS ■ 105 Ibs/acre 


RMS/standard deviation of YO ■ .3^ 
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that raportad In tha lltaratura. A curvlttnaar ralatlonship batwaan molstura 
avanab;Ilty at floMartng and taad yfald workad slightly batter than a simple 
linear one. It was, hOMavar, surprising that tha ralatlonship was not more 
significant. Certainly, tha buffer from stress provided by stored soil 
molstura during tha flowering period could have helped obscure tha YD to FP 
relationship. A soil molstura budget modeling approach, utilizing tha 
precipitation variables In conjunction with tha FCC (field capacity coefficient) 
variable did not produce significant results, unfortunately. 

The of the overall model was 0.90. If only environmental variables 
were used this value fell to 0.28. Thus for tha years of I98O and 1981, 
environmental variation over this geographic region apparently had little 
to do with sunflower yield variation. Vfhen Insect effects ware Included, tha 
root mean square was only 105 . This represents a 66 percent lnq>rovamant over 
the ability of the mean and standard deviation to characterize yields. Since 
yields ranged from 600 to 1800 Ibs/acre, the overall regression model represents 
a reasonbaly good fit. 

a 

On a discouraging note, the Importance of the 1C variable signifies a 
critical need to monitor or model insect populations associated with sunflower 
production. In the case of the seed weevil jnd sunflower midge, field 
monitoring is tedious and difficult, while modeling their phenology and 
degree of Infestation has not even been attempted to date. In remote sensing 
research, there Is some hope that the damage caused by insects through 
defoliation or premature senescence (early death of leaf tissue) might be 
detectable by multipectral sensors. However, damage caused by seed weevils 
or midge occurs in tha seed head and is not likely to be detectable by such 
sensors. In addition, bird damage occurs in a similar fashion and can be 
significant. Thus, it is important that some methods be developed to estimate 



34 


th« Impacts of those biological variables In order to batter model sunflower 
stress in certain environments. 

Because of the limited number of years (2) used in the CPM data base, 
the neture and degree of Impact from the environmental variables remains open 
to question. An analysis of yields from more diverse climates, either based on 
geographic distribution or long time series would be useful. The difficulty in 
this task Is finding data of sufficient detail to compute meaningful predictor 
variables. Future research should Include monitoring both growth stages and 
pests across a range of environments. Perhaps a cooperative effort between 
universities or counties interested in remote sensing applications to this 
crop would be more fruitful. 

Regression Analysis on the County Data Set 

Tables VII and VIII summarize the statistics from the best fit model for 
county sunflower yields. The most Important predictor variable, and the only 
one significant at P ■ .01 was CAPF. This confirms once again the importance 
of available moisture during the critical flowering period. This variable 
alone accounted for nearly 28 percent of historica'. yield variation. Over the 
time period studied, the CAPF variable ranged from 3.32 inches to only 0.39 
inches. In years with inadequate stored soil moisture, it would appear that 
a shortage of rainfall during the flowering period could sionificantiy reduce 
yields. According to the regression coefficient, each inch of rainfall during 
this period accounted for a nearly 200 Ibs/acre Increase In seed yield. 

The next important variable was TTV, as expected. Nearly ail agronomic 
crops show a positive time trend due to Improved genetics and cultural 
practices. The TTV variable accounted for 22 percent of the variation in yield, 
amounting to approximately a 23 Ibs/acre increase each year. It is difficult 
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T«blt VII 


A sutnmary of simple statistics from the best regression model for sunflower 
yield based on the county data base. 


Variable 

(units) 

Mean 

Value 

Standard 

Deviation 

Coefficient 
of Variation 

Maximum 

value 

Minimum 

Value 

CAYD (Ibs/acre) 

12A1 

265 

.214 

1780 

820 

TTV (coded) 

5.0 

2.62 

.524 

9 

1 

CAPO (Julian day) 

1A2 

7.41 

.052 

152 

131 

CAPF (Inches) 

1.73 

.67 

.388 

3.32 

0.39 

CATSF (days) 

4.83 

3.32 

.686 

16 

0 


A summary 
the county 

of regression statist 
data base. 

Table VI 1 
Ics for the 

1 1 

best sunflower yield 

model from 

Variable 

Regress ion 
Coefficient 

t-value 

Contribution 
to r2 

Stepwise 

r 2 

Intercept 

2223.85 

2.96 

- 

- 

TTV 

29.45 

2.44* 

.22 

.22 

CAPO 

- 9.57 

2.04* 

.065 

.285 

CAPF 

192.23 

4.14** 

.275 

.56 

CATSF 

-21.86 

2.14* 

.06 

.62 


Significant at P ■ .05 OF ■ 31 

** Significant at P • .01 
RMS - 173 


RMS/standard deviation of CAYD ■ .65 
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to know how significant • time trend verleble would be for other major sunflower 
producing areas. Based on this and other analyses, significant positive time 
trends In sunflower production have likely occurred In most major production 

areas. 

The CAPO and CATSF variables were approximately equal as predictors of 
sunflower yield. Both were significant at P ■ .0$. The planting date effects 
show an advantage to early planting, a concept discussed by other researchers 
In the literature. There was a 20 day range In county average planting dates 
for the 1971 to 1980 period. This compares with a 38 day range In planting 
dates from the CPH data base, where It was not a significant variable In the 
regression analysis. This fact may be due to the Importance of the Insect 
variable In the CPM analysis, which may have obscured any planting date effects. 

The CATSF was really somewhat of a pseudo-PET variable In the county 
analysis. The number of days during flowering when the dally temperature was 
90^F or higher may be Indicative of high potential evapotranspi ration (PET). 
Little is mentioned In the sunflower research literature about high temperature 
stress. The CATSF variable Indicated moderate reductions In seed yield with 
high temperatures during flowering. High teiripe natures were relatively absent 
In the CPH data set and were not used as a predictor. 

The for the county yield model was only 0.62, with an RMS of 173. 

This represents only a 35 percent Improvement over the mean and standard 
deviation in characterizing the sunflowrr yields of these counties. The 
relatively poor fit of this model may be due to the lack of more accurate 
growth stage estimates and to the absence of variables to account for pest 
damage. As sunflower culture has become more established In the Minnesota- 
North Dakota region, pests, likewise, have become more prevalent and had 
a greater impact on yields. 
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Th« GDD system, base degree F, may be a suitable method to estimate 
the phenology of sunfloMers. GOD estimates checked reasonably well with CPH 
surveys in cases where the hybrid maturity rating was known. Evaluating the 
GOO estimates of growth stage for the county data set was difficult due to 
the inadequacy of the USOA field surveys in the counties and to the range of 
maturity types represented in such data. As mentioned previously, the range 
of maturity in sunflower hybrids is rather narrow compared to field corn. 
However plant breeders have been working with this crop for a relatively short 
period of time. Therefore, the range in maturity types may Increase In the 
future as breeders make greater use of the genetic diversity in this crop. 

In this context. It is important to relate GOO estimates of phenology to 
specific sunflower maturity types. 

SUMMARY AND RECOHMENDAT I ONS 

These initial, albeit relatively simple analyses of both detailed field 
data and historical county data have provided some important evidence and 
guidelines for future modeling research on sunflowers. Some of the more 
important characteristics of sunflower production related to modeling stress 
include the following: 

1. There are significant time trends in sunflower yields since 1971. 

2. Similar to other agronomic crops, sunflowers appea.' to be most sensitive 

to temperature and moisture stress during the critical flowering (anthesis) 
and seed development stages. 

3. Using time trend, planting date, and environmental variables oriented to 

2 

the flowering stage as predictors of yields, an R of only 0.62 was 
obtained for Minn:*sota county data. 

4. Using detailed field survey data, including growth stage and insect records. 



38 


2 

•long with phenology*orlent«d •nvlronmental variables, an R of .90 was 
obtained on the CPh data base, where yields ranged from 600 to 1800 ibs/acre. 
5« As sunflower culture becomes more established, pest variables, especially 
insects, will likely become important predictors of yield. For the CPH 
study, the sunflower seed weevil and midge were very important variables. 

The morphological characteristics and nature of flowering in the sunflower 
crop would appear to invite remote sensing studies to determine the detectability 
of both growth stages and crop stress using muitispectral scanners. Such a 
study was initiated in Minnesota during the 1982 crop season. Remote sensing 
may have the potential to provide stress-related Information where environmental 
or biological variables affect the morphology, ontogeny, or color (necrosis or 
chlorosis) of the crop, as with drought, nutrient deficiencies or disease. 
However, some pests, such as the seed weevil and midge, may produce yield 
reducing effects not manifested by changes which are detectable by remote 
sensing techniques. For this reason, sunflower stress models should consider 
the use of meteorologically driven submodels for pest infestations. 
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Figure 2. Geographical Distribution of Complementary Aerial Photography 
Acquired to Assist In Landsat Analysis. 
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b) Two Seed Varieties in Same Field 


Normal Field 


c) Drainage Problem 


d) Soil Moisture Variability 


Figure 3* Sample Aerial Photographs Depicting "Normal" and Anomalous 
Field Cond i t i ons . 

* Original photographs are contained in the Final Report sent to Agricultural 
Applications Branch-SH3, NASA LBJ Space Center. 
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c) Moisture/Cultivator Damage d) Replant Condition 


Figure Sample Aerial Photographs Depicting Normal and Anamolous 
Field Cond i t i ons . 

* 

Original photographs are contained in the Final Report sent to Agricultural 
Applications Branch-SH3, NASA LBJ Space Center. 
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Crop water use of some crops commonly grown 
in the Northern Great Plains. 
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APPENDIX I 


0RIG!?4A1- FA'-' - ' -i 
OF rCC" QJAUW 


the data set listed below was derived from the university of MINNESOTA CROP 
PEST management FIELD SURVEYS FOR 1980 AND 1981, ALONG WITH DAILY WEATHER 
RECORDS KEPT BY NATIONAL WEATHER SERVICE COOPERATIVE OBSERVERS. THERE ARE 
A TOTAL OF 101 ROWS OF DATA, EACH WITH 16 DATA FIELDS, SEPARATED BY BLANK 
COLUMNS. the data FIELDS ARE DEFINED AS FOLLOWS: 


DATA 

FIELD 

1 

DATA 

FIELD 


DATA 

FIELD 

3 

DATA 

FIELD 

4 

DATA 

FIELD 

5 

DATA 

FIELD 

6 

DATA 

FIELD 

7 

DATA 

FIELD 

8 

DATA 

FIELD 

9 

DATA 

FIELD 

10 

DATA 

FIELD 

1 1 

DATA 

FIELD 

1£ 

DATA 

FIELD 

13 

DATA 

FIELD 

14 

DATA 

FIELD 

15 

DATA 

FIELD 

16 


FIELD IDENTIFIER (MISC. CODE) 

SAMPLE YIELD (LBS/ACRE) 

YEAR COLLECTED 

FIELD CAPACITY COEFFICIENT (SOIL TEXTURE) 
PLANTING DATE (JULIAN DAY) 

PRESEASON PRECIPITATION (INCHES) 

DATE OF ANTHESIS (JULIAN DAY) 

RAiNFALL 10 DAYS BEFORE TO 10 DAYS AFTER 
DATE OF ANTHESIS (STAGE A. 1) IN INCHES 
PET ESTIMATE FOR 5 DAYS BEFORE TO 5 DAYS 
AFTER BUD STAGE (STAGE 3.1) IN INCHES 
RAINFALL FOR 5 DAYS BEFORE TO 5 DAYS AFTER 
BUD STAGE ^ STAGE 3.1) IN INCHES 
PET ESTIMATE FOR 5 DAYS BEFORE TO 5 DAYS 
AFTER ANTHESIS (STAGE A. 1 ) IN INCHES 
RAINFALL FOR 5 DAYS BEFORE TO 5 DAYS AFTER 
ANTHESIS (STAGE A. 1) IN INCHES 
PET ESTIMATE FOR 5 DAYS BEFORE TO 5 DAYS 
AFTER SEED DEVELOPMENT (STAGE 5.1) IN INCHES 
RAINFALL FOR 5 DAYS BEFORE TO 5 DAYS AFTER 
SEED DEVELOPMENT (STAGE 5.1) IN INCHES 
FINAL PLANT POPULATION (IN THOUSANDS) 

INSECT CODE FOR SEED WEEVIL AND MIDGE 

0 = NO SIGNIFICANT DAMAGE 

1 = COUNTS APPROACHING ECONOMIC THRESHOLD 
£ = COUNTS EXCEEDING ECONOMIC THRESHOLD 


NASAl 101 15 FORMAT 


(bX, G4 

. 0, G3. 0, 

65. 2, 

G4. 0, G5. 2, G4. 

0, 7g: 

S. 2, G3. 0, G2. 0) 






10501 

1800 

80 

0. 80 

154 

5. 42 

218 

^ . .j3 

2. 74 

0. 34 

2. 18 

2. 42 

1. 77 

0. 20 

14 

0 

1£001 

1700 

80 

0. 80 

129 

7. 06 

207 

2. 55 

2. 23 

0. 77 

2. 27 

0. 37 

1. 80 

1.61 

15 

0 

13501 

1 650 

80 

0. 80 

129 

5. 42 

210 

4. 39 

2. 13 

0. 19 

2. 30 

1. 16 

1. 90 

0.28 

18 

0 

14001 

1200 

80 

0. 60 

136 

4. 58 

206 

3. 84 

2. 22 

0. 30 

2. 27 

0. 65 

1. 79 

0. 98 

17 

1 

1400£ 

1200 

80 

0. 60 

136 

4.58 

206 

3. 84 

2. 22 

0. 30 

2.27 

0. 65 

1. 95 

1. 35 

17 

1 

13003 

1300 

80 

0. 60 

134 

4. 58 

204 

3. 94 

2. 22 

0. 30 

2. 26 

1. 49 

1. 99 

0. 82 

21 

1 

15001 

0750 

80 

0. 75 

125 

8.98 

200 

2. 55 

2. 37 

1.03 

2. 62 

1. 34 

2. 39 

0.23 

19 

2 

16001 

1450 

80 

0. 60 

148 

4. 73 

216 

2. 47 

2. 61 

1. 70 

2. 62 

0. 26 

2. 15 

0. 39 

18 

0 

17001 

1800 

80 

0. 60 

128 

8.98 

206 

2. 18 

2. 60 

0. 79 

2.27 

0. 06 

1. 80 

0. 99 

15 

0 

18001 

1750 

80 

0. 75 

1 35 

8. 98 

206 

2. 18 

2. 60 

0. 79 

2. £7 

0. 06 

1. 80 

0. 99 

18 

0 

16002 

1600 

80 

0. 75 

135 

8. 98 

209 

1. 82 

2. 60 

0. 79 

2. 39 

0. 23 

1. 95 

1. 24 

20 

0 

19002 

1700 

80 

0. 75 

1 35 

4. 58 

209 

2. 57 

2. £9 

0. 09 

2. 37 

0. 65 

2. 01 

2. 50 

18 

0 

21001 

1650 

80 

0. 80 

133 

4. 10 

202 

2. 20 

2. 21 

0. 53 

2. 30 

0. 35 

1. 89 

1 . 

14 

0 

25001 

1650 

80 

0. 75 

135 

7. 06 

209 

2. 53 

2. 60 

0. 19 

2. 39 

0. 66 

1. 89 

0. 96 

1 7 

1^1 

28001 

1200 

80 

0. f 5 

136 

5. 76 

208 

2. 68 

2. 83 

0. 55 

2. 59 

0. 33 

2. 15 

1. 17 

14 

1 

30001 

120C 

80 

0. 75 

160 

4. 58 

204 

3. 94 

2. 28 

0. 30 

2. 26 

0. 63 

1. 92 

1. 35 

17 

1 

31001 

1200 

( g, 

0. 75 

159 

4. 58 

202 

3. 27 

2. 18 

0. 50 

2. 30 

0. 96 

1. 95 

1. 55 

15 

1 
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3 1 12103' 

1700 

80 

0. 75 

liv6 

4. 69 

198 

4. 05 

3. 17 

0. 31 


3. 36 

3.08 

0. 98 

14 

0 

330(211 

1800 

80 

0. 80 

135 

7. 06 

310 

3.47 

3. 60 

0. 19 

2. 34 

0.67 

1.89 

0.96 

18 

0 

33501 

1750 

80 

0. 60 

140 

5. 43 

309 

4.39 

3. 30 

0. 19 

3. 35 

1. 16 

1. '90 

1. 59 

19 

0 

34001 

1800 

80 

0. 80 

135 

7. 06 

303 

3. 55 

2.35 

0. 77 

2. 30 

0.05 

1.83 

1.49 

18 

0 

34003 

1750 

80 

0. 80 

135 

7. 06 

303 

3. 55 

3. 35 

0. 77 

3. 30 

0. 05 

1. 83 

1. 49 

16 

0 

35001 

1700 

80 

0.60 

130 

7.06 

305 

3.55 

3.36 

0. 77 

3.33 

0. 39 

1. 87 

1.61 

23 

0 

3&001 

1800 

80 

0. 65 

130 

5. 43 

310 

4. 39 

3. 18 

0. 19 

3. 34 

1.16 

1. 83 

0. 30 


0 

37001 

0670 

80 

0. 75 

158 

4. 00 

310 

6.97 

3. 17 

0. 30 

3. 30 

0. 89 

1.90 

1. 70 

15 


38001 

1750 

80 

0. 75 

143 

5. 43 

31 1 

4. 39 

3. 69 

0. ‘*3 

3. 33 

1. 30 

^ ^ Z’ 

X • «jc 

0. 37 

14 

0 

38003 

1700 

80 

0. 75 

143 

5. 43 

311 

4.39 

3. 69 

0. 43 

3. 33 

1. 30 

1 . .Jib 

0. 37 

1 6 

0 

39001 

1750 

80 

0. 80 

137 

4. 73 

303 

3. 19 

3. 51 

0. 51 

2. 47 

0. 91 

1. 84 

3. 1 1 

15 

0 

40001 

1800 

80 

0. 60 

135 

8. 98 

303 

3. 17 

3.36 

0. 73 

2.36 

0. 74 

1 . 83 

1.34 

1 ; 

(l 

40003 

1800 

80 

0. 60 

133 

7. 15 

301 

1.64 

3. 08 

3. 02 

3. 33 

1.01 

1. 73 

3. 43 

30 

0 

40003 

1800 

80 

0. 75 

135 

7. 06 

303 

3. 55 

3. 08 

0.08 

3. 36 

0.05 

1. 73 

1. 55 

30 

0 

43001 

1700 

80 

0. 60 

134 

4. 73 

310 

3. 33 


0. 40 

3. 34 

0. 16 

1. 80 

1. 37 


0 

46001 

1300 

80 

0. 80 

138 

4. 13 

304 

3. 66 

3.71 

1. 83 

3. 3& 

1. 02 

1.99 

3. 86 

33 

▲ 

47001 

1100 

80 

0. 75 

138 

4. 10 

306 

3.63 

3. 17 

0. 35 

3.37 

0. 67 

3. 05 

0. 92 

10 

1 

48003 

1400 

80 

0. 75 

137 

4. 73 

305 

3. 14 

3.89 

0. 77 

2. 46 

0. 83 

3. 33 

0. 70 

14 

0 

48003 

1500 

80 

0. 75 

139 

4. 73 

303 

3. 19 

3. 51 

0. 51 

2. 47 

0. 91 

3. w3 

0. 93 

14 

0 

51001 

0600 

80 

0. 75 

140 

4. 13 

309 

5. 63 

3.61 

0. 39 

3. 35 

0. 89 

1.85 

1. 70 

14 

3 

51003 

0630 

80 

0. 75 

143 

4. 13 

314 

6. 06 

3. 6l 

0. 39 

3. 38 

3. 87 

1. 85 

1. 70 

14 


53001 

0700 

80 

0. 75 

118 

4. 10 

199 

3. 44 

3.31 

0.53 

3. 48 

0.68 

2. 02 

0. 92 

16 

2 

5350 1 

1700 

80 

0. 75 

135 

8. 98 

303 

3. 17 

3. 36 

0. 80 

2. 36 

0. 74 

1.95 

1 *'4 

16 

0 

53503 

1 300 

80 

0. 75 

135 

8.98 

305 

3. 18 

3. 36 

0. 80 

r* *» *7 

^ ^ w 

0.52 

1.73 

... 4 

17 

1 

53503 

1400 

80 

0. 75 

135 

8. 98 

303 

3. 17 

3. 31 

0. 80 

3. 30 

0. 70 

. 

A . 

'4 

16 

1 

54001 

1 100 

80 

0. 65 

139 

4. 10 

313 

3.33 

3. 40 

0. 73 

3. 31 

1.07 

1 . ' 

.39 

14 

1 

55001 

1350 

80 

0. 75 
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5. 76 
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3. 44 

0. 63 

1 iZ 

1.61 

1^. ^ < 

0. 84 

18 

1 

56001 
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80 

0. 75 
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4. 13 

307 
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3. 38 

0.35 

3.39 

0. 84 

3. 1 1 

0. 99 

16 

0 

58001 

1700 
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0. 60 
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4. 73 

308 

3. 40 

3. 19 

0. 40 
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0. 16 
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1. 87 

15 

0 

59001 
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80 
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131 

4. 58 

304 

3.94 

2. 38 

0. 30 

3. 36 

0.63 

1.95 

1. 35 

30 

0 
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1500 

80 
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131 
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2. 17 

3. 40 
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3. 41 

0. 70 

1. 89 

1. 19 

17 

0 

60001 
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4. 00 
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3. 34 

0. 37 

3. 39 

0.67 

1.97 

3. 64 

30 

0 

60501 

1650 

80 
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3. 34 

0. 39 
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0. 83 

16 

0 

6300 1 
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3.30 

3. 17 
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3.30 
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3.05 

0. 93 

18 

0 

63001 

1600 

80 

0. 75 
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4. 73 

304 

3. 1 1 
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0. 77 

3. 49 

0. 91 
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1. 87 

18 

0 

63003 

1300 

80 

0. 75 

133 

4. 73 

301 

1. 88 

2. 51 

0. 51 

2. 59 

0. 91 

1.84 

3. 11 

13 

1 
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0. 03 
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0 
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0.51 
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0. 60 

133 

4. 73 
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PPPENDIX II 


THE YIELD DATP SET LISTED BELOW WPS DERIVED FROM USDP - 3RS YIELD DPTP FOR 
MPRSHPLl, polk, NORMPN, PND TRPVERSE COUNTIES OF MINNESOTA, COVERING THE 
PERIOD FROM 1971 THROUGH 1900. THE YEAR 1974 WPS ELIMINATED DUE TO 
SEVERE DROUGHT PND INACCURATE YIELD ESTIMATES FROM ABANDONED OR 
UNHARVESTED ACREAGE. THE PREDICTOR VARIABLES WERE DERIVED FROM USDA - 
CROP REPORTING SERVICE SURVEYS AND FROM LOCAL NATIONAL WEATHER SERVICE 
COOPERATIVE OBSERVER RECORDS OF DAILY TEMPERATURE AND PRECIPITATION. 

THERE ARE P TOTAL OF 36 ROWS OF DATA, EACH WITH 6 DPTP FIELDS, SEPARATED 
BY BLANK COLUMNS. THE DATA FIELDS ARE DEFINED AS FOLLOWS; 


DATA FIELD 1 

DATA FIELD 3 
DPTP FIELD 3 
DATA FIELD 4 


DPTP FIELD 


DPTP FIELD 6 


COUNTY CODE, WHERE 1 -MARSHALL, 2-POLK. 

3-NORMAN, AND 4-TRAVERSE 

YIELD IN LBS/ACRE (COUNTY AVERAGE) 

TIME TREND, WHERE 1=1971 AND 9-1980 
PLANTING DATE WHEN 30 PERCENT OF THE 
SUNFLOWER ACREAGE WAS SOWN ACCORDING 
TO USDA SURVEYS (JULIAN DAY) 

PRECIPITATION F.^OM 10 DAYS BEFORE TO 
10 DAYS AFTER THE ESTIMATED DATE OF 
ANTHESIS (USING A GDD METHOD) IN INCHES 
NUMBER OF DAYS WHEN THE DAILY MAXIMUM 
TEMPERATURES EXCEEDED 90 DEGREES F DURING 
THE PERIOD FROM 10 DAYS BEFORE TO 10 DAYS 
AFTER THE ESTIMATED DATE OF ANTHESIS 


COUNTY 

36 

6 FORMAT 
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